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In the present work, we have studied some features of the generalized Brans-Dicke (BD) model in which the scalar field is allowed 
to couple nonminimally with the matter sector. Extended holographic Ricci dark energy (EHRDE) has been considered in the 
above framework of BD cosmology. Some restrictions have been derived for the BD parameter co> and a stronger matter- chameleon 
coupling has been observed with the expansion of the universe. In this framework, the equation of the state parameter of EHRDE 
has behaved like quintom. Also, we have reconstructed the potential and coupling function for BD model for the EHRDE. It has 
been observed that the potential function is increasing as the matter- chameleon coupling is getting stronger. 



1. Introduction 

Accumulating the observational data of Supernovae Type la 
(SN la) by the year 1998, Riess et al. [1] in the High-redshift 
Supernova Search Team and Perlmutter et al. [2] in the Super- 
nova Cosmology Project Team have independently reported 
that the present universe is accelerating. The source for this 
late-time acceleration was dubbed "dark energy" (DE), which 
is distinguished from ordinary matter species such as baryons 
and radiation, in the sense that it has a negative pressure. 
Reviews on DE are available in [3-9]. The SN la observations 
have shown that about 70% of the present energy of the 
universe consists of DE [7]. The simplest candidate for DE is 
the so-called cosmological constant A, whose energy density 
remains constant [10] . If the origin of DE is not the cosmolog- 
ical constant, then one may seek for some alternative models 
for DE to explain the cosmic acceleration today. There are two 
approaches to construct models of DE other than cosmo- 
logical constant [7]: (i) to modify the right-hand side of the 
Einstein equations by considering specific forms of energy- 
momentum tensor T^ v with negative pressure; (ii) to modify 
the left-hand side of the Einstein equations. The representa- 
tive models that belong to the first category are the so-called 



quintessence [11], k-essence [12], and perfect fluid models 
[13] . The representative models that belong to the second class 
are "modified gravity" [14-22] models that include f(R) grav- 
ity [19], scalar-tensor theories [20], and braneworld models 
[21] . Among all scalar- tensor theories of gravity, the simplest 
one is the so-called Brans-Dicke theory (BD) introduced by 
Brans and Dicke [23] that modifies general relativity in accor- 
dance with Mach's principle. Cosmological models of classical 
BD theory were first studied in [24, 25]. In the subsequent 
decades, various aspects of BD cosmology have been widely 
investigated [26-34]. 

The present paper is aimed at working in chameleon 
Brans-Dicke cosmology. Thus, we discuss a bit of chameleon 
cosmology in this place. Khoury and Weltman [35, 36] 
employed self-interactions of the scalar field to avoid the 
bounds on such a field and dubbed such scalars to be "chame- 
leon fields" due to the way in which the field's mass depends 
on the density of matter in the local environment. In regions 
of high density, the chameleon "blends" with its environment 
and becomes essentially invisible to search for the violation of 
equivalence principle and the fifth force [35]. Reference [37] 
has shown that chameleon scalar field can provide explicit 
realizations of a quintessence model where the quintessence 
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scalar field couples directly to baryons and dark matter with 
gravitational strength. Reference [38] has illustrated that 
interacting chameleon field plays an important role in late- 
time universe acceleration and phantom crossing. Reference 
[39] proposed an interacting holographic dark energy model 
in chameleon-tachyon cosmology by interaction between the 
components of the dark sectors. Based on two independent 
functions of the scalar field, [40] constructed an exact 
solution describing the evolution of the type Bang-to-Rip 
with the phantom divide line crossing in the chameleon 
cosmology. 

Now, we come to the chameleon Brans-Dicke cosmology. 
In the models, where a nonminimal coupling between the 
scalar field and matter system is considered by introducing 
an arbitrary functions of the scalar field, the scalar field is 
regarded as as a chameleon field [33]. Reference [41] was the 
first to introduce this kind of chameleon-matter coupling in 
the BD model to achieve accelerated the expansion of the 
universe. In a recent paper, [33] considered a generalized BD 
model with power-law form of the scale factor and the cou- 
pling functions as the inputs and allowed nonminimal cou- 
pling with the matter sector to show that accelerated expan- 
sion of the universe can be realized for a contrained range of 
exponents of the potential function. Treating the BD scalar 
field as a chameleon scalar field and taking a nonminimal 
coupling of the scalar field with matter, [42] studied cos- 
mological implication of holographic dark energy in the BD 
gravity. Reference [43] investigated the cosmological appli- 
cations of interacting holographic dark energy in BD theory 
with chameleon scalar field which is nonminimally coupled 
to the matter field and observed that phantom crossing can be 
constructed if the model parameters are chosen suitably. For 
different epochs of the cosmic evolution, [44] investigated the 
BD chameleon theory of gravity and obtained exact solutions 
of the scale factor and scalar field and potential. Reference 
[45] discussed late-time dynamics of a chameleonic general- 
ized BD cosmology with the power law chameleonic function 
/(0) oc 0", where n is a real parameter motivated from string 
theories. 

Making an approach to dark energy in the frame of 
quantum gravity is the well-known holographic dark energy 
(HDE) inspired by the holographic principle [46]. The HDE, 
whose density is p A = 3c 2 M 2 L~ 2 was proposed by Li [47]. 
Subsequent studies on HDE from various points of view 
include [48-51]. In the present work, we are considering a 
special form of HDE [52] dubbed as "extended holographic 
Ricci dark energy" (EHRDE) [53] , whose density has the form 



p A = 3M 2 p (aH 2 + pH), 



(1) 



where M 2 is the reduced Planck mass and a and /3 are con- 
stants to be determined. Like [42], we will assume that p A and 
pressureless dark matter are conserved separately, and we will 
assume a nonminimal coupling between the scalar field and 
the matter field. 



2. EHRDE in Chameleon BD Cosmology 

We begin with the chameleon BD theory in which the scalar 
field is coupled nonminimally to the matter field via the 
action [33] 

S = l - I d 4 x V=^ (j>R - ^V,, W - 2V + 2/ (4>) L m ) . 

(2) 

In (2), R and denote the Ricci scalar field and BD scalar field, 
respectively. The /(</>) and V(<fi) are analytic functions of the 
scalar field. The matter Lagrangian density, denoted by L m , 
is coupled with via the function / "(0). This function allows 
a nonminimal coupling between the matter system and the 
scalar field. If /(0) = 1 , we get back the BD action with poten- 
tial function [33]. Varying the action with respect to the 
metric g^ v and 0, one gets the field equations [33] 

0G, v = ^ v + /(0)T™ (3) 

(2o> + 3) U(j> + 2 (2V - V'<f) = T m f - 2/> m , (4) 

where □ = V^V^.V^ represents the covariant derivative, T m = 
g hv T™, and the prime denotes differentiation with respect to 
0. In (3), 



+ (V fJ V^-^ v D0)-V(0)^ v) 
-2 <HV=#LJ 



T 



(5) 



(6) 



F9 <V V 

Because of the explicit coupling between matter system and 0, 
the stress tensor T™ v is not divergence free. 

Now, we will apply the above framework to a homoge- 
neous and isotropic universe described by the Friedman- 
Robertson-Walker metric 

ds 2 = -dt 2 + a (t) ( -^—r + r 2 dQ 2 \ . (7) 



-1,+1 or 



l-kr 2 

The universe is open, closed, or flat according to k 
0. In a spatially flat universe, (3) and (6) yield 

2 f cv 2 V 

F 2 2 



(8) 



3(H + H 2 ) = 



3p 



0(2o) + 3) 

; 2 



y0f' + (a>(y+^) + 1 )/ 



1 

2 2co + 3 



3V + (2cv - 3) 



V 



(2co + 3) (0 + 3H0) - 2 (2V - 0V') 
= p((l-3y)/ + 2y0/'). 



(9) 



(10) 
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In (8) and (10), p = p m + p A and y = p A /(p m + p A ) (since dark 
matter is pressureless, p m = 0). This approach is similar to 
[42, 43]. 

In the present work, we will investigate two cases. In one 
case, we will assume special forms for V", /, a, and and 
derive conditions that strengthen matter- chameleon coupling 
with the expansion of the universe. In another case, we will 
not assume any form for V and /. Rather we will reconstruct 
them for the EHRDE. 

3. Discussion 

3.1. Case I. In this section, we will consider a set of ansatz 
for the potential function, analytic function, scalar field, and 
scale factor. Based on them, we will determine constraint on 
the BD parameter co. Subsequently, based on the constraints 
on co and other parameters, we will investigate the behavior of 
the EHRDE in BD chameleon cosmology. Following [33], we 
choose 



a (t) = a t n ; 



/(</>) = fo^ 

<l>(t) = <l> t m . 



(11) 

(12) 



Since M 2 = 1/8ttG and in BD theory <j> oc G 1 [43], we 



can take EHRDE in the chameleon BD as 
p A = 3c/> (ccH 2 + fiH) . 



(13) 



And using conservation equation p A + 3H(p A + p A ) = 0, we 
get 

p A = ^(l - 3H 2 - 2H) a + (H - H (h + 3HH))^ . 

(14) 



Using (11)-(14) in (8), we have 



,mU-m+2-3n , ,2-m+mL , ,mU 



(15) 



where q v £ 2 , £ 3 , and £ 4 are terms involving constant terms 
of ansatz (11). In (15), we observe that the right-hand side is 
constant with t involved in the left-hand side. This is possible 
if the powers of t are 0. In (15), ml 2 ± 0. However, its coefficient 
3nf o 0Q (na - /3)(3 + co + 3y(/ 2 + co)) can be set equal to 0. Since 
na-fi = will lead to p A = 0, we can set 3 + co+3y(l 2 + co) = 0. 
Finally, we get the restrictions 



3(i + y/ 2 ) 

CO = , 

1 + 3y 
-2 = m (l 2 - l) - 3n, 
-2 = m(l l -l). 



(16) 



(17) 



The above restrictions are also found valid if we use the ansatz 
(11) and (13)-(14) in (9). In (16), we have for the accelerated 
expansion of the universe 3y + 1 > 0. Hence, (16) finally leads 
to a constraint on the BD parameter as 



U 



co > 



1 + 3y 



(18) 




Figure 1: Evolution of EoS parameter y against t. The solid, dot- 
dashed and dashed lines have the parameter combinations a = 
2.3, 1.3, 1.5, p = 4.5, 2.5, 1.5, l Y = 0.3, 0.5, 0.2, and l 2 = 4.5, 6, 7. 



In 1973, co > 5 was consistent with the known data. By 1981, 
co > 30 was consistent with the known data [54] . In 2003, evi- 
dence derived from the Cassini-Huygens experiment shows 
that the value of co must exceed 40,000 [55]. Hence, we can 
rewrite 



co > 



l + 3y 

Hence, U > 4. From (17), we have 



» 1. 



n = 



m = 



(19) 



(20) 



If \ x > 1, then m < 0, which implies that the scalar field is 
a decreasing function of t. Since n has to be positive, l x > l 2 
from (20). Again, if \ x < 1, then m > 0. Since l 2 > 0, ml 2 > 0. 
Thus, for l x < 1, coupling function / oc t mXl increases with 
time, that is, matter- chameleon coupling gets stronger as the 
universe expands. 

Choosing appropriate parameters (keeping the restric- 
tions in mind), we have plotted the EoS parameter y = 
pip for different combinations of parameters in Figure 1. 
The EoS parameter clearly shows a transition from y > -1 
(quintessence) to y < -1 (phantom). Hence, the EoS param- 
eter is found to exhibit "quintom"-like behavior. Next, we plot 
in Figure 2 the deceleration parameter q [53] 



1 3p 
- + — . 

2 2p 



(21) 



We observe that the deceleration parameter has a clear tran- 
sition from positive to negative side. This indicates the tran- 
sition from decelerated to accelerated phase of the universe. 
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Figure 2: Evolution of deceleration parameter q (21) against t. The 
solid, dot-dashed and dashed lines have the parameter combinations 
a = 2.3, 1.3, 1.5, p = 4.5, 2.5, 1.5, l x = 0.3, 0.5, 0.2, and l 2 = 4.5, 6, 7. 



mn (4 + m) (3 + 2co) t Hypergeometric2Fl 



2 + ra 6 (n - 1) + m 

~> r, 



_ 4 - 2m - 6n 2 (3n + m - 2) 
ftp (m ~ 2) 0o (na - /3) 

PmO 

+ (2 + m) (- 4n 2 - 2mn (2 + n) 

-m 2 (1 + 2n) + m 3 (1 + a))) 
x Hypergeometric2Fl 



4 + m 6n - 8 + m 



. 4 - 2m - 6n ' '2 (3n + m - 2) ' 



a 3 , (m - 2) O (na - /3) 

PmO 



(24) 



Using (24) in (22), we get the reconstructed potential 
function as 



3.2. Case II. In this section, we will only assume a(t) = a t n 
and 0(f) = (j>Qt m . We will use these ansatz to reconstruct the 
potential V and the analytic function / of the chameleon BD 
cosmology On putting these ansatz in (8), we get an equation 
that involves V, f and the remaining terms are functions of t 



T7- 1 .-In r 

V = -^Pmot f 



m-2 



(22) 



+ — [-m 2 o) + 6n (m + n - (n - /3) /)] . 



We again differentiate it with respect to t and get a new equa- 
tion involving V and /. Using this in (9) along with the said 
ansatz, we get the following differential equation on / with t 
as the independent variable: 

df (t) 3 (m + In) 



dt 



2t 



-fit) 



o 3 / ,m+3rc-l 

= 3a o o f 

x [mn {t (3 + 2co) - 4 - In) ( 23 ) 
+m 3 (1 + co) - m (1 + In) - An) 

x (2 (Pmot 2 ~ 4 (fit - 2) O (na - /3)))~\ 
Solving (23), we get reconstructed analytic function /(0) as 

/(0) = /(*bO 



= r 



m/2 



Pmo (8 + 6m + m 2 )f 2 



2a^ (2 + m) (4 + m) p m0 
x (- (2 + m) (4 + m) p m0 t 2 
x(2C lPmO ^ )/2 + a o 3 o 
x ( - 6mn + 6n 

x(-n + Clt 3/2(m+2n) (na-?)) 
+ m 2 a;) ) ) 
+ 4 (Pmo^ 2 + 3a 3 n0 o £ m+3n (na - /?)) 

mn (4 + m) (3 + 2co) t Hypergeometric2Fl 

2 + m 2 + m 

1, , 1 + , 

4 - 2m - 6n 4 - 2m - 6n 

a 3 (m-2)0 o (na-£)r 2+m+3n " 

PmO 

+ (2 + m) (-4n 2 - 2mn (2 + n) 

-m (2n + 1) + m (1 + a))) 
x Hypergeometric2I 7 l 

4 + m 



4 + m 

1, , 1 + 

4 - 2m - 6n 



4 - 2m - 6n 



a 3 , (m - 2) O (na - /3) f 

PmO 



-2+m+3n 



(25) 
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Figure 4: Reconstructed V (0) against for various combinations of 
parameters. 



In Figures 3, 4, and 5, we have plotted the reconstructed 
/(0) against 0, reconstructed V((f)) against 0, and recon- 
structed V(0) against /(0), respectively We observe that 
/(0) -> Oas0 -> Oand^(0) -> as -> 0. These indi- 
cate the satisfaction of one of the sufficient conditions for real- 
istic models. In Figure 5, we observe that the potential func- 
tion V((f>) is increasing with the coupling function /(0). Thus, 
we infer that as the matter-chameleon coupling gets stronger, 
the potential function also increases. In Figure 6, we have 
plotted /(0) against t for a range of positive values of m. This 




Figure 6: Reconstructed /(0) against t for a range of values of m. 



figure shows that /(0) in an increasing function of the cosmic 
time t. This indicates that with the expansion of the universe, 
the matter- chameleon coupling gets stronger, where the 
coupling function has been reconstructed for EHRDE. 

4. Conclusion 

In the present work, we have studied some features of the 
generalized BD model in which the scalar field is allowed to 
couple nonminimally with the matter sector. Extended holo- 
graphic Ricci dark energy has been considered in the above 
framework of BD cosmology. The energy density has been 
considered in the form p A = 30(aH 2 + to accommodate 
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the BD scalar field in the dark energy density The work has 
been carried out in two different perspectives. Firstly, we have 
assumed the ansatz V{$) = V </> Zl ; /(</>) = / o ?2 ; a(t) = a t n ; 
and (f)(t) = (f) t m and derived some restrictions on the BD 
parameter co. Under this restriction, we have investigated 
the equation of state parameter for the dark energy and 
a quintom-like behavior has been observed. Furthermore, 
the deceleration parameter has shown a transition from 
decelerated to accelerated phase of the universe. In this 
scenario, it has been observed that the restrictions are leading 
to a stronger matter- chameleon coupling with the expansion 
of the universe. Secondly, we have considered only a(t) = a t n 
and <fi(t) = (j) t m and reconstructed the coupling function and 
potential for the extended holographic Ricci dark energy in 
BD model. We have observed that the potential is increasing 
with the increase in the matter- chameleon coupling. Further- 
more, we have seen that for any positive value of m, /(</>) is an 
increasing function of t. This means that for this reconstruc- 
tion, the coupling function is increasing; that is, the matter- 
chameleon coupling is getting stronger with passage of cos- 
mic time. 
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